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Pure gravitational plane waves are considered as a special case of spacetimes with two commuting 
spacelike Killing vector fields. Starting with a midisuperspace that describes this kind of spacetimes, 
we introduce gauge-fixing and symmetry conditions that remove all non-physical degrees of freedom 
and ensure that the classical solutions are plane waves. In this way, we arrive at a reduced model with 
no constraints and whose only degrees of freedom are given by two fields. In a suitable coordinate 
system, the reduced Hamiltonian that generates the time evolution of this model turns out to vanish, 
so that all relevant information is contained in the symplectic structure. We calculate this symplectic 
structure and particularize our discussion to the case of linearly polarized plane waves. The reduced 
phase space can then be described by an infinite set of annihilation and creation like variables. We 
finally quantize the linearly polarized model by introducing a Fock representation for these variables. 
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I. INTRODUCTION 

The analysis of gravitational waves has received lately 
a great deal of attention, both from the theoretical and 
the experimental points of view B. In particular, there 
is an increasing interest in the study and quantization 
of spacetimes with two commuting Killing vectors that 
describe wave solutions in source-free general relativity 
[2-5]. One of the reasons for this interest, apart from 
possible applications in cosmology and astrophysics, is 
that this kind of spacetimes provides a good arena to test 
quantization techniques and discuss conceptual issues in 
gravity. The symmetry reduction of Einstein gravity in 
the presence of two commuting Killing vectors leads in 
general to models with an infinite number of degrees of 
freedom, i.e., to midisuperspace models. It is commonly 
accepted that the quantization of these models, which is 
given by a true quantum field theory, may be significant 
for discussing basic features of the outstanding theory 
of quantum gravity. In contrast, minisuperspace trunca- 
tions of general relativity, like, e.g., most of the gravi- 
tational systems quantized in the literature ||, are too 
simple to capture the field complexity that should be 
present in full quantum gravity. Another motivation for 
studying the reduction of general relativity by two com- 
muting spacelike Killing vectors is that it admits a similar 
formulation to that for coset space cr-models coupled to 
two-dimensional gravity and a dilaton j7|. The quanti- 
zation of any of these two types of systems should hence 
allow a better understanding of the quantum physics of 
the other. 

The spacetimes with two commuting Killing vectors 
that have been quantized so far are the family of cylin- 
drically symmetric gravitational waves and the Gowdy 
model with the spatial topology of a three-torus. The 



Gowdy universes are vacuum spacetimes with compact 
sections of constant time [§). They can be thought of 
as inhomogeneous spacetimes filled with gravitational 
waves. The quantization of the Gowdy model when the 
topology of the sections of constant time is that of a 
three-torus was analyzed in Ref. 0] . Preliminary studies 
of this quantization had been carried out by Berger || 
(assuming that the Killing vectors of the model were hy- 
persurface orthogonal), and by Husain and Smolin JlQ] . 
On the other hand, the quantization of cylindrical grav- 
itational waves with linear polarization was achieved by 
Ashtekar and Picrri ||, completing previous works on 
the subject by Kuchaf [Q and Allen Q. The effects 
produced by quantum fluctuations in this cylindrically 
symmetric model were studied in Ref. |l2[ ]. Finally, the 
more general case of cylindrical waves with two polariza- 
tions was recently quantized by Korotkin and Samtleben 

In addition to possessing two commuting spacelike 
Killing vector fields, the above families of spacetimes sat- 
isfy another condition: the existence of two-surfaces that 
are everywhere orthogonal to the group orbits spanned 
by the Killing vectors. For spacetimes of this class, a 
specially useful metric function is the square root of the 
determinant of the two-metric that corresponds to the 
group orbits |l3|]. We will denote this positive function 
by W. The spacelike, timelike, or null character of the 
gradient of W is invariant under coordinate transforma- 
tions. In this paper, we will study the quantization of 
another example of gravitational waves that belongs to 
this class of spacetimes, namely, the case of pure grav- 
itational plane waves [13|,|l4j. Our analysis will exhaust 
the discussion of the different posibilities for the gradi- 
ent of W, in the sense that this gradient is timclikc for 
the Gowdy model with the topology of a three-torus B, 
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spacelike for cylindrical gravitational waves |3j , and null 
for plane waves. 

Gravitational plane waves are a special type of pp- 
waves, i.e., plane- fronted gravitational waves with paral- 
lel rays |l3l , |l5| ]. They were first studied by Baldwin and 
Jeffrey |16| , and admit a five-dimensional group of isome- 
tries with a three-dimensional Abelian subgroup that acts 
on null hypersurfaces |T?j| . Gravitational plane waves 
can be interpreted as describing the gravitational field at 
great distances from finite radiating bodies |ll|,[l6| . Some 
of these waves are related to certain limits of Bianchi cos- 
mological models and inhomogeneous generalizations of 
such models fl4| , |l8|| . On the other hand, the collision 
of two gravitational plane waves has received intensive 
study, and a lot of solutions that represent this collision 
are known [ pTpg| ] . 

In this work, we will restrict our considerations to the 
case of pure gravitational plane waves, namely, plane 
waves that are solutions to Einstein equations in vac- 
uum. Preliminary studies of the quantization of space- 
times with the symmetry of these waves have been car- 
ried out by Neville |2^,^l| and Borissov p^ j, mainly in 
the context of the Ashtekar formalism for gravity. Here, 
we will start with the geometrodynamic formulation for 
general relativity, and discuss the quantization of pure 
gravitational plane waves as a particular case of space- 
times that possess two commuting spacelike Killing vec- 
tors. 

The rest of the paper is organized as follows. The geo- 
metrical properties of pure gravitational plane waves are 
reviewed in Sec. II. With a suitable choice of (local) co- 
ordinates, we cast the metric of this family of solutions 
in a convenient (3+1) form. This metric can be written 
in terms of two arbitrary functions, and corresponds to 
spacetimes with two commuting Killing vector fields that 
are subject to certain conditions. In Sec. Ill, we consider 
the Hamiltonian formulation for spacetimes that admit 
two commuting spacelike Killing vectors and perform a 
partial gauge fixing. The reduction of the system is com- 
pleted in Sec. IV, where we impose gauge-fixing and sym- 
metry conditions that remove all non-physical degrees of 
freedom and guarantee that the classical solutions are 
precisely pure gravitational plane waves. We determine 
the reduced phase space of the model obtained in this 
way and its symplectic form. It is also shown that, in 
an appropriate set of coordinates, the dynamical evolu- 
tion of this reduced model is trivial, so that all relevant 
information about the system is in fact contained in the 
symplectic structure. In Sec. V, we particularize our 
study to the family of linearly polarized plane waves in 
vacuum gravity. We prove that the symplectic form on 
the reduced phase space of this model can be regarded as 
that corresponding to an infinite set of annihilation and 
creation like variables. We then proceed to quantize the 
model by introducing a Fock representation. Finally, we 
summarize our results and conclude in Sec. VI. 



II. PURE GRAVITATIONAL PLANE WAVES 

Gravitational plane waves are solutions to Einstein 
field equations that possess a five-parameter group of 
isometries |17| ]. For these spacetimes, the metric can al- 
ways be written in the form UjA: 



ds 2 



-dUdV + H ab (U)X a X b dU 2 



[dX 



a \2 



(2.1) 



where we have employed the index notation a,b = 1,2. 
The coordinates {U^V^X 1 ^ 2 } are called harmonic co- 
ordinates, and run over the whole real axis. It is easy 
to check that dy corresponds to a c ovariantly constant 
null vector field, so that metrics (2.1) are a special class 
of pp- waves |13|. On the other hand, the only non- 
vanishing component of the Ricci tensor for these metrics 
is Ruu = — {H\i +i?22)- Solutions to Einstein equations 
in vacuum are hence given by all symmetric matrices H a b 
with vanishing trace. In this case, the plane wave is said 
to be purely gravitational. 

Although the use of harmonic coordinates for plane 
waves has clear advantages, because they allow to cover 
the totality of the spacetime with a single chart and ex- 
press the Ricci tensor in terms of a simple metric func- 
tion, the spacetime symmetries become much more ap- 
parent in the so-called group coordinates. In such coor- 
dinates, the metric reads [fD&Ol: 



ds 2 



-dUdV + h ab (U)dx a dx b 



(2.2) 



Here, V, x , and x 2 run over the entire real line, and 
the coordinate U has in general a restricted range, as we 
will see below. Obviously, d x i and d x 2 are two commut- 
ing spacelike Killing vector fields for the above metrics, 
and the group orbits spanned by them are everywhere 
orthogonal to the surface with coordinates U and V. It 
is convenient to write the two-metric h ab of these group 
orbits as 



h u = e z -y' 2 , h 12 = -ve z ^' 2 , 

h 22 = {v 2 + e y)e z - y ' 2 . 



(2.3) 



Hence, the condition that h ab is positive definite is auto- 
matically satisfied provided that the functions v, y, and z 
are real. Note that the determinant of h a b is equal to e 2z . 
Thus, from our discussion in the Introduction, we have 
W = e z and, since this function depends only on the co- 
ordinate U, its gradient is null with respect to the metric 
( pT^ ) , as we had anticipated. Finally, the functions v and 
y describe, respectively, the diagonal and non-diagonal 
degrees of freedom of the metric e~ z h ab , which has unit 
determinant. 

In group coordinates, the only non-trivial Einstein 
equation in vacuum is 



R 



uu 



I ab (Phab 

"2 dU 2 



1 dh ab dh, 



4 dU dU 



0. 



(2.4) 
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with h ab being the inverse of the two- metric h a b- Substi- 
tuting Eq. (pTsj) in this expression, we get 



In (z'o) 



d 2 (e z ' 2 ) 
dU 2 



e z/2 

"16" 



dy 
dU 



(2.5) 



Therefore, W 1 / 2 = e z / 2 must be a convex function of 
the coordinate U . Since this function is non-negative, it 
must then generally vanish at some point. The metric 
becomes degenerate at that point and a coordinate sin- 
gularity appears. As a consequence, the spacetime for 
pure gravitational plane waves cannot be described by a 
unique chart in group coordinates |T^| . 

Actually, the fact that W 1 ^ 2 is a non-negative convex 
function implies that it must vanish at either one or two 
points in the allowed (simply connected) interval of vari- 
ation for the coordinate U. The only exception is pure 
Minkowski spacetime, where v, y, and z are constant. 
When W x l 2 has two zeros, it increases from the first of 
them until a maximum is reached. This maximum can 
always be made equal to the unity by means of a suit- 
able rescaling of the real coordinates x a . On the other 
hand, when W 1 / 2 has a single zero, it can be chosen as 
a strictly increasing function of U by performing a re- 
versal of the null coordinates U and V, if necessary. In 
this case, W 1 / 2 will either tend to infinity or to a finite 
value in the limit of large positive U's. If the limit is 
finite, an appropriate rescaling of x a sets it equal to the 
unity. Hence, we see that in all cases there exists a region 
of the spacetime where the function W 1 ^ 2 increases with 
U and ranges in the interval (0,1). In this region, z is 
strictly increasing from minus infinity to zero. In addi- 
tion, although Minkowski spacetime has apparently been 
excluded from our discussion, it turns out that the above 
considerations apply as well to flat spacetime, because 
the solution to Eq. (2.5) in which v and y are constant 
and W 1 ! 2 is linear in U can be cast in Minkowskian form 
by a (local) coordinate transformation JL7[ . 

In the region where z increases in (— oo,0), we can 
introduce a change of coordinates 



dU = e 2 



^y' 2 du 



(2.6) 



so that the function z gets fixed in terms of our new 
coordinate u. For instance, we can make 



z(U[u}) = -e- u = z (u). 



(2.7) 



In Eq. fl2.6| ), w is an unknown function of u. This coor- 
dinate is assumed to run over the entire real axis. Our 
change of coordinates is then well-defined and invertible, 
because zq is strictly increasing for u G M and has the 
range Ft~ . With this change, the coordinate singularity 
at e 2 = has been driven to minus infinity. Notice also 
that our ignorance about the original function z(U) is 
enco ded in w. However, the vacuum Einstein equation 
( [2.5| ) becomes now a first-order differential equation in w 
that can be solved explicitly. The solution is 
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7 [0/) 2 + 4(t/) 



Wo, 



(2.8) 



where uo is any fixed real constant and the prime stands 
for derivative with respect to u. The constant of integra- 
tion that should appear in this formula has been absorbed 
by a scale transformation in the coordinate V. 

If we finally define a time via V — 2t — u, we arrive 
at the following line element for pure gravitational plane 
waves: 

ds 2 = e 2w °- Zo+v/2 [-dt 2 + {du - dt) 2 ] + h ab dx a dx b , (2.9) 



with h a b the metric obtained from Eq. (2.2) by setting 
z = z - All coordinates in this expression are rea l and 
have unrestricted ranges. Recalling relations ( |2.7| ) and 
( [2.8| ), we see that this family of vacuum solutions have 
only two degrees of freedom, which are given by the func- 
tions v and y of the coordinate u. 

These are the spacetimes that we will study in the 
rest of this paper. As we have shown, the above metric 
describes the most general pure gravitational plane wave 
solution, with the only caveat that it does not correspond 
to the whole of the spacetime that can be covered with 
harmonic coordinates. The coordinate system that we 
have adopted is in fact quite similar to that employed in 
Rcf. 1 17 1, the main difference being that we have chosen 
zo = —e~ u and allowed u to take all real values, while zo 
was made equal to 2 In |u| in Ref. [jlTj. In this latter case, 
u has a restricted range [e.g., u € (0, 1) for zq 6 M~]. 

We can regard Eq. (2.9) as providing the standard 3+1 
decomposition of the metric for spacetimes with two or- 
thogonal surfaces, but with the lapse and (non- vanishing 
component of) the shift vector particularized to take the 
values N = e ^-z /2+ y /i and N u = _ L In addition, 

the spacetimes considered admit two commuting space- 
like Killing vector fields, d x i and d x z, and satisfy another 
condition, namely, that the metric functions are indepen- 
dent not only of x a , but also of the time coordinate t. In 
the Hamiltonian formalism for general relativity, there is 
a simple way of imposing the time independence of the 
metrics (2.9) on classical solutions, once the existence of 
the two spacelikc Killing vectors has been assumed. We 
first note that zo is a fixed function of u, and that wo 
is determined by the vacuum Einstein equations to be 
given by Eq. 



(2. 



in terms of v and y. It hence suffices 
to ensure that the time derivatives v and y vanish on 
classical solutions. Employing then the definition of the 
extrinsic curvature p3| , with x l = {x 1 , :e 2 , u}, it is 
not difficult to check that, for the metrics under discus- 
sion, the requirements v = y = on classical solutions 
are equivalent to the relations 

Ah l/2 e -2z 0+V /2 Kii = 2z , _ y / ) 

Ah l/2 e -2z +y/2 Ki2 = vy , _ 2vz , o _ 2^ ( 2 

where h = e 2w o+^o+y/2 j g determinant of the induced 
metric hij. Moreover, these relations turn out to imply 
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the vanishing of v and y even when the values of Wq, the 
lapse function, and the shift component N u are changed 
in the metrics (2.9), provided that N u does not depend on 
the coordinates x a and the quotient N/N u = — /i 1 / 2 e _z ° 
is not modified. 



III. SPACETIMES WITH TWO KILLING 
VECTOR FIELDS 

In this section, we will discuss the Hamiltonian formu- 
lation for spacetimes that possess two commuting space- 
like Killing vectors, paying a special attention to the case 
of pure gravitational plane waves. 

We will only consider spacetimes that admit a global 
foliation, so that the metric can be written in the (3+1) 
form: 



ds 1 



N 2 dt 2 + h i:i (dx< + N'dt)^ + N J dt). (3.1) 



For the coordinates of the sections of constant time, we 
will employ the notation {x 1 } = {x l ,x 2 ,u\. In addition, 
we will assume that the coordinates t and u run over 
the real line, and impose that d x a (a = 1,2) are the 
two spacelike Killing vector fields. As a consequence, the 
metric must be independent of the coordinates x a . This 
implies that the area S = J dx 1 dx 2 appears as a constant 
global factor in the gravitational Hilbcrt-Einstein action 
and in the symplectic form for geometrodynamics. We 
can always absorb this area in Newton's constant. In 
order to do this, renormalization is required when the 
area S is infinite |^0| (this case can be considered as the 
limit S — ► oo of the compact case). For convenience, we 
further set the effective Newton's constant obtained in 
this way equal to l/8n. 

The fundamental Poisson brackets and momenta cano- 
nically conjugate to hij are then given by 



{h l3 (u),U kl (u)} = 6 { *S l )s(u- u), 
IF = ^h 1/2 (h ik h jl - h l3 h kl )K kl , 



(3.2) 
(3.3) 



where 8* is the Kronecker delta, S(u) is the Dirac delta on 
the real line, the indices in parenthesis are symmetrized, 
and is the inverse of the induced metric. Thus, the 
symplectic form on phase space is 



dudlT 3 Ad/ij 



(3.4) 



with the symbols d and A denoting the exterior derivative 
and product, respectively. 

The gravitational system has two kinds of constraints: 
the Hamiltonian constraint H. and the momentum con- 
straints Tli, which generate spatial diffeomorphisms. In- 
stead of dealing with H, we will consider its densitizcd 
version, H = h x l 2 H. In terms of phase-space variables, 
these constraints take the expressions 



H 



(3) 



R + {hikhji + h a h jk - h 2j h k i)U ij U kl = 



Hi = -2/i«D fc Ir^ = 0. 



(3.5) 



Here, Di is the covariant derivative compatible with the 
induced metric and ^ R its curvature scalar. The time 
derivative of any function on phase space is then 



./ = d t f + {/, J du (NH + N'Hi)}, 



(3.6) 



where dt is the partial derivative with respect to the ex- 
plicit time dependence and N = hT x l 2 N is the densitized 
lapse function. 

Using that the metric does not depend on the coordi- 
nates x a , it is not difficult to check that the momentum 
constraints H. a can be rewritten as 



Ha 



-2 h a ,n l 



(3.7) 



For the globally hyperbolic spacetimes with two commut- 
ing spacelike Killing vectors that we are studying, it is 
then possible to fix the gauge freedom associated with 
these constraints by demanding that, on the sections of 
constant time, the w-line is orthogonal to the group or- 
bits spanned by d x a . This requirement can be rephrased 
in the form 



0. 



1,2. 



(3.8) 



To see that the above conditions lead to a well-posed 
gauge fixing, we must prove that, in general, the gauge 
orbits generated by the constraints H. a intersect trans- 
versely t he s urface that is defined on phase space by con- 
ditions (3.8) and the constraints of the system. This is 
equivalent to prove that the gauge-fixing conditions have 
non- vanishing Poisson brackets with H. a . A straightfor- 
ward calculation shows 



{K 



,JduN b H b } = h ab (N b y. (3.9) 



Taking into account that h ao is positive definite, we con- 
clude that these brackets are generally different from zero 
(it suffices to take (N a )' ^ 0). 

We must also check that our gauge fixing is compati- 
ble with dynamics, in the sense that there exists a choice 
for the components iV a of the shift vector (which are the 
Lagrange multipliers of 7i a ) so that, on the constraint 



surface, the gauge conditions fl3.8j ) are preserved by the 
dynamical evolution. In other words, modulo constraints 
and gauge-fixing conditions, h au must vanish for a par- 
ticular choice of iV a . Recalling the independence of the 
metric on x a , one can see that 



h ab (N b Y + 4Nh ab l 



rb u 



(3.10) 



where use of the gauge conditions has been made after 
calculating Poisson brackets. On the other hand, the 
solution to H a — is given by 
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IF" = h ab f b (t), 



(3-11) 



with f a being two arbitrary functions of t that must be 
real, because IT™ G JR. Thus, our gauge fixing is in fact 
preserved by the evolution provided that 



(N a y = -mh uu h ab f b (t). 



(3.12) 



Let us now integrate both terms in the above equality 
over it, and contract the result with f a {t). If we assume 
that the components N a of the shift vanish in the limit 
of large w's, we get 



/ duNh uu h ab f a f b = 0. 



(3.13) 



It is worth commenting that, together with the condi- 
tions h au = 0, the requirement liniu^-i-oo N a — can 
be viewed as ensuring that the group orbits spanned by 
d x a are asymptotically orthogonal to the surfaces with 
coordinates t and u. This requirement should hence be 
satisfied at least in spacetimes where those surfaces are 



everywhere orthogonal. Returning to Eq. (3.13), we note 



that, since the integrand in that expression is real, it has 
to vanish at least at one point if the integral is equal to 
zero. Taking into account that the two-metric Nh uu h ab 
is positive definite, it follows that f a = at that point. 
The functions f a must then be identically zero, because 
they are independent of u. Therefore, one obtains that, 
when N a vanishes for large absolute values of u, the so- 
lutions to the constraints TL a — and consistency condi- 
tions h au = are 



IP" = 0, N a = 0. 



(3.14) 



The vanishing of h au and 7V a implies that the space- 
time metric is block-diagonal, so that the assumption 
linin^-too A^ a = employed in our gauge fixing turns 
out to guarantee that the surface with coordinates x a is 
actually orthogonal to that described by t and u. 

After this partial gauge fixing, the canonically con- 
jugate pairs (h aU jTl au ) and the momentum constraints 
TC a are eliminated from the system. The remaining geo- 
metrodynamic momenta and components of the induced 
metric provide a canonical set of variables on the phase 
space of the gauge-fixed model. The line element of this 
reduced model has the form 



ds 2 = e 2w - z+ yl 2 [-e 2z N 2 dt 2 + (du + N u dt) 2 ] + 

e z - y l 2 [(dx x ) 2 - 2vdx 1 dx 2 + {v 2 + e y ){dx 2 ) 2 ]. (3.16) 

Here, all metric functions are allowed to depend on the 
coordinates u and t, and we have densitized the lapse 
using that 



2w+z+y/2 



(3.17) 



Similarly to the situation found for plane waves, the con- 
dition that the induced metric is positive definite is now 
easily imposed: it suffices to demand that the functions 
q a are real. On the other hand, since the change of vari- 
ables performed is just a point transformation, it is pos- 
sible to find a set of momenta {p a } that are canonically 
conjugate to our new metric variables. These momenta 
are 



_ yruu dk uu 

PQ dq<* 



n 



ab d h ab 



dq° 



(3.18) 



where h uu a nd h nb are regarded as functions of q a . Us- 
ing relation (3.3), we then obtain that, in terms of the 
extrinsic curvature and our metric variables, 

Pv =-h 1 ' 2 e-*-y' 2 (vK 11 +K 12 ), 

p w = -h^ 2 e- z - y / 2 [(v 2 + e y )K xl + 2vK 12 + K 22 ], 
1, 



= -h}/ 2 e- 2w+z -yl 2 K n 



(3.19) 



Finally, the gravitational model considered is still sub- 
ject to two constraints, which are given by the restriction 
of 7i u and 7i to the case in which the metric is indepen- 
dent of the coordinates x a and h au — n au = 0. After 
some trivial calculations, one can show that the momen- 
tum constraint TC U takes the form 



74 = -Pw + Pw W ' + PvV' + PyU' + P Z Z', 



(3.20) 



and that the densitized Hamiltonian constraint has the 
expression 

W = 44 Kv'f + 4{v') 2 e-y - Az'{Aw' + y'- bz') + 16z"] 
16 



p 2 v e v - 2p w p y - p w p z + Ap z y . 



(3.21) 



ds 2 = -N 2 dt 2 + h uu {du + N u dt) 2 + h ab dx a dx b , (3.15) 

where the metric functions depend only on the coordi- 
nates t and u. 

At this stage, it is convenient to carry out a change of 
metric variables from h uu and h a b to a set of functions 
{q a } = {v, w, y, z} that is analogous to that employed to 
describe the family of metrics (2.9). In this way, our dis- 
cussion will be straightforwardly applicable to the analy- 
sis of pure gravitational plane waves. The change is such 
that the line element adopts the expression 



IV. GAUGE FIXING AND SYMMETRY 
REDUCTION 

We will now apply the analysis of the previous sec- 
tion to the study of pure gravitational plane waves. We 
will complete the gauge fixing of the system and intro- 
duce symmetry conditions that ensure that all classical 
solutions are plane waves in vacuum, so that the reduced 
mod el th at is obtained really describes the family of met- 
rics (pj|). 
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Let us first remove the gauge freedom associated with 
the constraint H Ul which generates diffeomorphisms in 
the coordinate u. A way to fix this gauge is to provide 
a one-to-one correspondence between the coordinate u 
and a metric variable. This can be done in the case of 
the plane wave metrics that we are considering, because 
the variable z can then be chosen as a strictly increasing 
function of w Si?, namely, the function zq defined in 
Eq. (2.7). In order to fix the gauge freedom, we then 
introduce the condition 



Xu 



0. 



(4.1) 



Since z' ^ 0, the constraint H u — can be easily solved 
to find the momentum canonically conjugate to z: 



Pz=p° z = -r(p' w 



p w w' - p v v' - p v y'). (4.2) 



Therefore, if our gauge fixing is admissible, we can elim- 
inate the canonical pai r (z ,p z ) from the system. 

Actually, condition (11) leads to a well-posed gauge 
fixing, because (modulo that condition) 



{ Xu , / duN u H u } 



z' o m 



(4.3) 



Obviously, this Poisson bracket differs from zero when 
N u 7^ 0, because z' does not vanish at any point. On 
the other hand, if our gauge condition is compatible with 
the dynamical evolution, there must exist a choice for ./V" 
(i.e., the Lagrange multiplier of H u ) such that 

Xu = {Xu,J du {NH + N U H U )} 

= - Pw N + z' N u = 0. (4.4) 

In the second line of this equation, we have substituted 
z — zq. The condition \u = is thus preserved by the 
evolution if 



N u = 



PwN 



(4.5) 



Note that the component N u of the shift given by this 
formula is always well-defined, because zq never vanishes. 
We hence conclude that the gauge fixing introduced is 
acceptable. It is worth remarking that, to arrive at this 
result, the only property of the function zq that has been 
used is that its derivative with respect to u is everywhere 
finite and different from zero. In this sense, the explicit 
form of this function, as well as its range, are irrelevant. 

After this gauge fixing, we obtain a reduced model 
whose phase space is described by the canonical set of 
variables {v,w,y,p v ,p w ,p y }. The symplectic form is 

SI = J du (dp v A dv + dp w A dw + dp y A dy), (4.6) 

and the system is subject to the densitized Hamiltonian 
constraint 



H r — H(z = z 0} p z — p z ) = 0. 



(4.7) 



In order to complete the reduction and attain a midis- 
uperspace model for pure gravitational plane waves, we 
still need to introduce two types of conditions. Firstly, 
we have to impose that all classical solutions of the sys- 
tem are independent of the time coordinate, so that these 
s olu tions correspond to the family of plane wave metrics 
(pTgj) . Secondly, we have to find a condition that results 
in removing the gauge freedom associated with the densi- 
tized Hamiltonian constraint. Both tasks can be achieved 
in the following way. 

We hav e show n at the end of Sec. II that, for metrics of 
the form ( |3.16| ) with z — z and N u = —e z °N, the time 
independence of the metric functions v an d y on classi- 
cal so lutions is ensured by relations ( 2.10 ). Using Eqs. 
( 3 . 1S| ) , those relations can be expressed as the conditions 



Xi = Py ~ 
X 2 = 2e»->„ 



-(y'-2z , )=0, 
= 0. 



(4.6 



Taking into account identity ( |4.5| ), we see that these con- 
ditions guarantee that v and y do not depend classically 
on time provided that 



Xo = p w + z' Q e z " = 0. 



(4.9) 



When this is the case, it is clear that conditions (4 



will also imply the time independence of the momenta p v 
and p y . In addition, the momentum p w will not depend 
on time if the requirement xo = is satisfied. The only 
phase-space variable that might then be time dependent 
on classical solutions is w. This possible time dependence 
can nevertheless be eliminated by demanding that 



X 3 = w - wo = 0, 



(4.10) 



where u>o is the function defined in Eq. (2.S). In this 
way, the function w is fixed to take the same expression 
in terms of v and y as it adopts in pure gravitational 
plane waves. On the other hand, the constraint of the 
model can be rewritten after some manipulations as 

H r = (2y' X i + e*°-V X 2 + 2u/ X o - 2 X ' )^ + 4 X ? 

zz 

- 2(xi + e zo X3)Xo + e zo X 'o- (4-11) 

So, once the symmetry conditions xi — X2 = arc 
imposed, the requirement xo = is just a solution to 
H r = 0. Therefore, the surface defined on phase space 
by xi — {I = 0, ...,3) is simply a section of the con- 
straint surface. In the following, we will show that the 
reduction to this section is consistent and that, in the 
reduction process, the condition X3 = allows us to fix 
the gauge freedom associated with H r - We will employ 
the symbol w for weak identities, namely, identities that 
are satisfied modulo xi — 0- Besides, we will denote the 
set {H r , xi, X2, X3> by {x^}, with A = 0, 3. 







To prove that the reduction is admissible, we first have 
to show that the conditions and constraint \A — consti- 
tute a second-class system on the section of the constraint 
surface that we are analyzing. Let us define 



x(g) 



J A 



(4.12) 



where {g A } is any set of functions of the coordi- 
nate u, i.e., functions of u that are infinitely differentiable 
and have compact support. Clearly, imposing xa = is 
equivalent to demand that x(g) vanishes for all choices 
of the functions g A . The set {xa} 1S then second-class 
provided that no combination of the form x(d) (other 
than the zero constant) commutes weakly under Poisson 
brackets with all of the xa' s - Using expression (4.11) for 
the constraint, and taking into account that all terms in 
that expression are quadratic in xi except the last one, 
we obtain 



{n r ,x(g)} 



-e z °g' 3 . 



(4.13) 



Since e z ° ^ and the function g 3 is Cfi°, it follows that 
x(g) commutes weakly with TL r if and only if g% vanishes. 
Let now x(g) be f°r S3 — 0- After some calculations, 
one gets 

2 1 

E 9a{xa, x(g)} « - g (e zo gi + \^ y - z »g\)'- (4.14) 

A=\ 

Given that the term in parenthesis is a sum of non- 
negative functions and that gi , 52 £ Cq° , the vanishing 
of the Poisson brackets of x{g) with xi and xi turns out 
to imply that g\ and g 2 must be the zero function. On 
the other hand, it is not difficult to check that 



{X3,x(ff)}«-(e Zo 3o)' 



(4.15) 



where x(g) denotes the restriction of x(g) to the case in 
which go is the only function in {gj.} that differs from 
zero. Recalling that g G Cq°, we see that the above 
Poisson bracket vanishes weakly provided that go = 0. 
Therefore, x(ff) must be exactly equal to zero in order 
to commute weakly with the set {x^}, as we wanted to 
prove. 

Let us now show that the section of the constraint sur- 
face determined by xi =0 (/ = 0, 3) is preserved by 
the dynamical evolution for a su itabl e choice of the den- 
sitized lapse. From expression (4.11) and the fact that 
the variables {p a } = {v, y,p v ,p w ,p y } commute with xo, 
we get 



{p a , / duNH r } w 0. 



(4.16) 



As a consequence, the variables p a turn out to be time 
independent on classical solutions. To deduce this result, 
one actually needs not demand that X3 vanishes. More- 
over, one can then straightforwardly see that the sym- 
metry conditions xi = X2 = are preserved in time, and 



that, when these conditions are imposed, the solution to 
the constraint given by xo = is invariant under the 
evolution, regardless of the value taken by the densitized 
lapse. On the other hand, employing again the expres- 
sion for the constraint and that {x3,Xo} = { w ,Xo}, we 
obtain 



X3 « -{e z °N)' 



w. 



(4.17) 



So, compatibility of the condition X3 — with the dy- 
namical evolution ensures that the variable w does not 
depend classically on time, and implies 



N = F{t)e 



(4.18) 



The arbitrary function F that appears in this expression 
can in fact be absorbed by means of the time redefini- 
tion t = J F(t )dt. In this way, we arrive at N = e~ z ° . 
From Eqs. ( |4. 5| ) and (|L9|), the only non- vanishing com- 
ponent of the shift vector is then given by N u = — 1. We 
hence see that the reduction proposed for the system is 
acceptable and that the condition X3 = allows us to fix 
the value of the densitized lapse function. Therefore, this 
condition removes the gauge freedom associated with the 
constraint TL r . In addition, we note that, in agreement 
with our previous discussion, all phase-space variables 
are indeed independent of the time coordinate on classi- 
cal solutions after reducing the system to the section of 
the constraint surface defined by xi = 0. 

In this reduction, the pair (w,p w ) is eliminated via 
the gauge-fixing condition X3 = and the solution to 
the constraint xo — 0, while the symmetry conditions 
Xi = Xi = determine the momenta p v and p y in terms 
of the variables v and y. The reduced model that is 
attained is free of constraints and possesses only two de- 
grees of freedom, which are described by the metric vari- 
ables v and y. The line ele ment for this reduced model 
can be obtained from Eq. (3.16) by substituting z = zq, 
w = wq , and the values obtained with our gauge fixing for 
N u and the densitized lapse function. This lin e clement 
coincides formally with that given in Eq. (2.9) for pure 
gravitational plane waves. The only difference is that the 
functions v and y of the coordinate u may now depend 
also on time. However, this time dependence is ruled 
out on classical solutions, as we have shown above. The 
metric may depend on time only when non-classical tra- 
jectories are allowed, e.g., in a quantum theory. Besides, 
it is clear that the set of classical solutions for our reduced 
model is precisely the family of gravitational plane waves 
in vacuum considered in Sec. II. As a consequence, our 
model describes in fact that family of plane waves. 

Obviously, the dynamics of this reduced model is triv- 
ial, because the variables v and y remain constant in time 
under the classical evolution. The reduced Hamiltonian, 
which generates the time evolution via Poisson brackets, 
must hence vanish. Of course, the same conclusion is 
reached by considering the reduced action for the system. 
This action can be obtained from the Hilbert-Einstein ac- 
tion in the following manner. We adopt the (3+1) decom- 
position of the metric explained in Sec. Ill, impose that 
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the metric is independent of x a (a — 1, 2), and absorb the 
area that corresponds to these coordinates in Newton's 
constant, setting the resulting effective constant equal to 
1/87T. Since all gravitational constraints are eliminated 
in the reduction process, we can make them equal to zero 
in the Hamiltonian form of the action. The reduced ac- 
tion S r that we are seeking can then be computed from 
the Lagrangian density IP^/iy. Substituting H au = 0, 
relations (3.18), and the values taken in our model by 



the variables [w, z,p v ,p w ,p y ,p z ), we get 



dtdu —(y'y + Av'e v v). 



(4.19) 



In arriving at this result, we have disregarded surface 
terms that are evaluated on sections of constant time, 
for they do not modify the reduced Hamiltonian (nor 
affect the symplectic structure). We see that the reduced 
action is linear in time derivatives. This implies that the 
reduced model has indeed a vanishing Hamiltonian in the 
set of coordinates adopted. 

Given that the dynamics is trivial, all relevant informa- 
tion for quantizing the system is encoded in its symplectic 
structure. The symplectic form f2 r on the phase space of 
the re duc ed model is provided by the pull-back of the 
form (4.6) to the surface determined by the conditions 
Xi = (I = 0, ...3). Taking into account that xo = 
ensures that p w is a fixed function of the coordinate u 
(because so is z ), we obtain 



Sir 



du— [dy' A dy + 4d(w'e" 2 ') Adv\. (4.20) 



Actually, this for m co uld have also been deduced from the 
reduced action ( ff,19| ). The quantization of the system 
can then be attained by quantizing the algebra of Pois- 
son brackets that follows from this symplectic form. In 
principle, the algebra seems rather complicated. This sit- 
uation contrasts with that found in other midisuperspace 
models which describe spacetimes with two commuting 
spacelike Killing vectors, like, e.g., the Gowdy model 
where the Poisson algebra obtained after (an almost com- 
plete) gauge fixing is simple, although the complexity of 
the system shows up in the expression for the reduced 
Hamiltonian. Nonetheless, it turns out that there ex- 
ists at least a class of pure gravitational plane waves for 
which this Poisson algebra becomes manageable, namely, 
the case of linearly polarized plane waves. From now on, 
we will restrict our discussion to this type of waves. 



V. LINEARLY POLARIZED PLANE WAVES 

Pure gravitational plane waves for which the metric 
function v vanishes are called linearly polarized ]ljj] . One 
can actually impose the vanishing of v as an additional 
symmetry on the midisuperspace model for plane waves 
in source-free gravity that we have studied. This can be 



done, for instance, by including the condition Xi = v = 
in the process of gauge fixing and symmetry reduction 
explained in the previous section. Similar arguments 
to those presented for general plane waves in vacuum 
gravity show that the set formed by Xi an d XA (with 
A = 0,...,3) is second-class on the section of the con- 
straint surface defined by the requirements xi = Xi = 
(I = 0, ..3). Here, we have employed the notation intro- 
duced in Sec. IV. On the other hand, we have already 
seen that the metric variable v is invariant under the 
dynamical evolution when the conditions xi — ar e sat- 
isfied, so that the symmetry condition Xi — is compati- 
ble with dynamics. Therefore, the reduction obtained by 
demanding that xi and Xi vanish is fully consistent. 

The reduced model that one obtains in this manner 
is such that its phase space can be described by a single 
metric variable: the field y. As it happens for plane waves 
with two polarizations, this variable depends on the co- 
ordinate u and, although it might in principle evolve in 
time, its time dependence is forbidden on classical so- 
lutions. It is easy to check that the action S r and the 
symplectic form O r for this red uced model coin cide in 
fact with those given in Eqs. (LIE) and (4.20) when 
particularized to the case v = 0. Defining 



Y = 



c.zo/2 



u 



2V2 



we can then write 



S r = / dtduY'Y, n r = / dudY' A dy, 



(5.1) 



(5.2) 



where we have used that zq is a fixed function of u and 
neglected surface contributions to the action that come 
from sections of constant time. Note that the reduced ac- 
tion S r corresponds indeed to a system with zero Hamil- 
tonian and symplectic structure given by n r . 

Like in the case of waves with two polarizations, the dy- 
namics of the reduced model for linearly polarized waves 
is trivial in the set of coordinates adopted. The function 

Y is constant in time and the reduced Hamiltonian van- 
ishes. Hence, in order to quantize the system, we only 
need to quantize the algebra of Poisson brackets that fol- 
lows from the symplectic form f2 r . 

On the other hand, employing that the coordinate u 
runs over the whole real axis, we can express the function 

Y as the Fourier transform 



Y 



1 



dk 



(a k e +a k e ). 



(5.3) 



The complex functions a k and a* kl with k 6 M + , provide 
then a complete set of variables on the phase space of 
the reduced model. These functions of k might depend 
as well on time; however, since Y is time independent on 
classical solutions, they turn out to be classical constants 
of motion. Substituting the above expression for Y in the 
symplectic form, we obtain 
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dk dat A dak- 



(5.4) 



As a consequence, the only non- vanishing Poisson brack- 
ets between the phase-space variables ak and a* k are 



{^,4} = -iS(k-k). 



(5.5) 



In addition, the fact that the function Y is real implies 
that a\ must be the complex conjugate of ak- Therefore, 
we can interpret the algebra of Poisson brackets for our 
model as corresponding to an infinite set of harmonic 
oscillators, described by the annihilation and creation like 
variables ak and a£, respectively. 

The quantization of this algebra can be performed by 
standard methods. For instance, one can simply intro- 
duce a Fock representation |25|]. In that case, our phase- 
space variables are represented as annihilation and cre- 
ation operators, ak and a£, and the Hilbert space of phys- 
ical states can be constructed by the repeated action of 
the creation operators on a vacuum, which is destroyed 
by all the operators ak- The algebra of commutators is 
given by the quantum analog of Eq. (5.5), 



(5.6) 



where we have set % = 1, and k, k 6 M + . Assuming that 
the vacuum has unit norm, the inner product is totally 
determined by requiring that the relations under complex 
conjugation between the variables ak and a* k are realized 
in the quantum theory as adjointness relations between 
operators. In particular, it is then possible to find nor- 
malized states which are formally similar to the states 
with n particles or to the coherent states of ordinary 
quantum field theory in flat spacetime. 

In this way, one attains a consistent and well-defined 
mathematical framework for analizying the quantum 
physics of spacetimes with two commuting Killing vec- 
tors that describe pure gravitational plane waves with 
linear polarization. 



VI. CONCLUSIONS AND FURTHER 
COMMENTS 

Pure gravitational plane waves have been considered as 
a special class of spacetimes with two commuting space- 
like Killing vector fields. We have discussed the structure 
of the reduced phase space and the quantization of this 
family of solutions to vacuum general relativity. 

We have first seen that, with a suitable choice of coor- 
dinates, the metric for these plane waves can be cast in a 
3+1 form in which the values of the lapse and shift are 
fixed and all metric functions depend just on a coordi- 
nate of the sections of constant time, namely, the coordi- 
nate u E M. The other two coordinates of these sections 
(which run over the real axis) correspond to Killing vec- 
tors, with group orbits that are orthogonal to the surface 
described by u and the time coordinate. We have shown 



that the determinant of the metric for these group or- 
bits can be chosen as a strictly increasing function of 
the coordinate u, with range given by the interval (0, 1). 
In addition, we have solved explicitly Einstein equations 
in vacuum, and introduced metric variables whose real- 
ity ensures that the induced metric is positive definite. 
The family of metrics that we have obtained in this man- 
ner represent the most general plane wave solution in 
source-free gravity, with the only caveat that they do 
not describe the whole spacetime which can in principle 
be covered with harmonic coordinates. These metrics are 
determined by two arbitrary functions of the coordinate 
u, i.e., the functions v and y. 

With the aim at studying this family of metrics, we 
have considered the Hamiltonian formalism for space- 
times that admit two commuting spacelike Killing vec- 
tors. Assuming the spatial topology to be that of M 3 , 
it has been proved that the condition that there exists a 
surface orthogonal to the group orbits removes the gauge 
freedom related with diffeomorphims of the coordinates 
which correspond to the Killing vectors. We have then 
particularized our analysis to the case of pure gravita- 
tional plane waves. Using that the determinant of the 
metric for the group orbits is an increasing function of the 
coordinate u, we have been able to eliminate the degrees 
of freedom associated with diffeomorphisms of the u-line. 
In order to complete the gauge fixing and guarantee that 
the classical solutions of the system are plane waves, we 
have next introduced two types of requirements. On 
the one hand, we have imposed the symmetry conditions 
(I.J), which can be interpreted as relations between the 
metric and the extrinsic curvature that are satisfied for 
plane waves. Once these conditions arc included, the 
densitized Hamiltonian constraint, which is the only con- 
straint that still remains on the system, is easily solved. 
Moreover, all metric variables are then classically time 
independent, except possibly the variable w, defined by 
means of Eq. (3.16). On the other hand, we have de- 
manded that w adopts the same expression in terms of 
v and y as it does in pure gravitational plane waves. We 
have shown that, in this way, the gauge freedom associ- 
ated with the densitized Hamiltonian constraint is totally 
removed and the time dependence of w is ruled out on 
classical solutions. 

After this process of gauge fixing and symmetry reduc- 
tion, we have arrived at a midisupcrspace model which 
is free of constraints and whose line element coincides 
formally with that for plane waves in vacuum general 
relativity, although the two functions v and y that deter- 
mine the metric may depend not just on the coordinate 
u, but also on time if non-classical trajectories are al- 
lowed. The phase space of this reduced model has only 
two degrees of freedom, which are described by the vari- 
ables v and y. These variables remain constant in time 
on all classical solutions, so that the dynamical evolu- 
tion on phase space is given by the identity map. This 
implies that the reduced Hamiltonian vanishes in the set 
of coordinates employed. The same conclusion has been 
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reached by computing the reduced Hilbert-Einstein ac- 
tion. Up to surface contributions on sections of constant 
time, this action is linear in time derivatives, so that it 
indeed corresponds to a system with vanishing Hamilto- 
nian. Since the dynamics is trivial, all the information 
needed for quantizing the reduced model is provided by 
its symplectic structure. Starting with the symplectic 
form for spacetimes with two spacelike Killing vectors, 
and calculating the pull-back to the section of the con- 
straint surface determined by the gauge-fixing and sym- 
metry conditions imposed on our model, we have then 
obtained the symplectic form on the reduced phase space. 

There exists at least a case in which the algebra of 
Poisson brackets that follows from this symplectic form is 
simple enough as to allow a straightforward quantization, 
namely, the case of linearly polarized plane waves. The 
restriction to this subfamily of plane waves in source- 
free gravity has been achieved by demanding that the 
function v vanishes. This requirement can be viewed as 
an additional symmetry condition on the system. The 
phase space of the resulting reduced model is described 
by the metric variable y or, equivalently, by the set of 
variables {cifc,a£; k € M + }. These variables are classical 
constants of motion obtained from the Fourier transform 
of the product of y and a fixed function of the coordinate 
u. We have finally proved that the algebra of Poisson 
brackets for the reduced model of linearly polarized plane 
waves can actually be understood as corresponding to an 
infinite set of harmonic oscillators whose annihilation and 
creation like variables are and a* k , respectively. The 
quantization of the model can then be readily performed 
by introducing a Fock representation for these variables. 

The mathematical framework constructed in this man- 
ner can be used to study the quantum physics of pure 
gravitational plane waves with linear polarization. In 
particular, one can try to define operators for the met- 
ric components and analyze the quantum fluctuations of 
the geometry. Another appealing possibility is to dis- 
cuss the role played in the quantum theory by coherent 
states and see whether these are peaked around classical 
solutions, comparing the results with those obtained for 
linearly polarized waves with cylindrical symmetry p|JT^] . 
The consideration of these issues will be the subject of a 
future research. 

Another direction for further investigation consists in 
introducing a scalar field in the model and studying its in- 
teraction with gravitational plane waves that are linearly 
polarized. This interaction has already been discussed 
in the general case of waves with two polarizations using 
a semiclassical approximation, i.e., considering only the 
quantization of the scalar field |^4j . It would be interest- 
ing to investigate whether the conclusions reached with 
this semiclassical analysis are valid to some extent in a 
purely quantum theory. 

Let us conclude with a couple of remarks about our re- 
duced model for linearly polarized plane waves. The first 
comment refers to the algebra of Poisson brackets that 
we have found. It is not difficult to see that, at a given 



instant of time, the values of the metric variable y at two 
different points do not generally commute under Poisson 
brackets. This implies that, in the quantum theory, there 
cannot exist states in which the variable y takes a well- 
defined value at every point on a section of constant time, 
even though the points on this section are spacelike sep- 
arated. In fact, this result is not so surprising. We have 
seen that the reduced Hamiltonian of the system vanishes 
in the set of coordinates adopted, so that the quantum 
evolution is dictated by the identity operator. Therefore, 
if it were possible to determine completely the value of 
y on & section of constant time, the same would happen 
for the entire spacetime, including points with timelike 
separation. The second of our remarks concerns the van- 
ishing of the reduced Hamiltonian. This conclusion is 
valid in the set of coordinates that we have employed to 
describe the spacetime. In another coordinate system, 
however, the Hamiltonian of our reduced model can gen- 
erally differ from zero. Suppose, for instance, that we 
change coordinates from u to x = u — t. Usi ng the ex- 
pression for the reduced action given by Eq. (p.2|) in the 
linearly polarized case, one can prove that the reduced 
Hamiltonian of the system becomes then 



dx(d x Y) 2 = / dkkakOf., 



(6.1) 



where Y is the variab le d efined in Eq. ( |5.1[ ) , and we have 
substituted relation ( |5.3| ) after replacing the coordinate 
u with x. In the new coordinate system, the interpre- 
tation of our reduced model as a collection of harmonic 
oscillators applies hence not only to the Poisson algebra, 
but also to the Hamiltonian that generates the time evo- 
lution. 
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